Stochastic gravitational wave background from cold dark matter halos 
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The current knowledge of cosmological structure formation suggests that Cold Dark Matter 
(CDM) halos possess a non-spherical density profile, implying that cosmic structures can be po- 
tential sources of gravitational waves via power transfer from scalar perturbations to tensor metric 
modes in the non- linear regime. By means of a previously developed mathematical formalism and a 
triaxial collapse model, we numerically estimate the stochastic gravitational-wave background gen- 
erated by CDM halos during the fully non-linear stage of their evolution. Our results suggest that 
the energy density associated with this background is comparable to that produced by primordial 



tensor modes at frequencies u ~ W — 10 Hz if the energy scale of inflation is V 



1/4 



1-2x10' 



GeV, and that these gravitational waves could give rise to several cosmological efi^ects, including 
secondary CMB anisotropy and polarization. 



PACS numbers; 98.80.Cq; DFPD 04/A-18 



I. INTRODUCTION 

Sources of gravitational waves (GW) are commonly 
separated in two types: astrophysical and cosmological. 

The first kind of sources can produce a stochastic back- 
ground which provides interesting information on the dis- 
tribution of compact objects at relatively low redshifts, 
such as star formation and supernova rates, black-hole 
growth mechanisms and other important phenomena. 
Such a background is generated by neutron stars, black 
holes and the associated binary systems, which emit in 
the frequency range Ve ~ 10^ — 10'*Hz. (e.g [1, 2]), or by 
galactic merging of unresolved binary white dwarfs with 
frequencies in the range Ve ~ 10~^ — 10~^Hz [3-7]. 

Besides binary systems of super-massive black holes 
in the galaxy center, which could emit at ~ 
hence detectable by LISA {e.g. Ref. [8]), the principal ex- 
ample of gravitational waves of cosmological origin is rep- 
resented by the relic radiation which has been generated 
by quantum fluctuations of the metric tensor during the 
inflationary era. The detection of this relic background 
would shed light on the physics of the very early Universe, 
since its strain amplitude is proportional to the square of 
the inflation energy scale. Primordial backgrounds can 
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be generated by various mechanisms and are character- 
ized by a large frequency interval which extends from a 
few 10^^* Hz to a few GHz, allowing their detection by 
markedly different ways of observation [1]. 

One of the best strategies for detecting the relic grav- 
itational radiation is to exploit the imprints it leaves on 
the Cosmic Microwave Background (CMB) temperature 
anisotropy and polarization [9-11]. More specifically, the 
CMB photons are very sensitive to primordial GWs with 
frequencies ~ 10~^''Hz, which correspond to the comov- 
ing size of the Hubble radius at last scattering, when ten- 
sor metric modes, being damped by the horizon entering, 
produce the largest amount of temperature quadrupole 
anisotropy and, consequently, by Thomson scattering, 
the largest amount of polarization [12]. It may be shown 
that the curl component in the polarization pattern, com- 
monly known as B-mode, is excited by vector and tensor 
cosmological perturbations only; therefore, if initial fluc- 
tuations are created very early, e.g. during inflation so 
that the vector growth is damped, primary B-modes can 
be produced only by tensor perturbations and, therefore, 
a possible detection will represent the incontrovertible 
proof of their existence [13-15]. 

Unfortunately, there are mechanisms that can produce 
secondary B-modes, the principal one being represented 
by the gravitational lensing, i.e. cosmic shear (CS) [16], 
which distorts the primary CMB pattern, in particular 
converting E- into B-modes [17]. Luckily, although com- 
parable, B-modes from primordial GW exhibit their peak 
at multipoles I « 100, corresponding to the degree scale. 
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while, for lenscd B-modcs, the peak is at I w 1000, cor- 
responding to the arcminute scale. Nonetheless, if the 
energy scale of inflation is V^/^ < 2 - 4 x 10i'''GeV, the 
CS-induced curl is a foreground for the / ~ 50 — 100 
primordial GW-induced B-polarization. This important 
contamination has to be removed in order to detect 
relic gravitational waves [18]. However, as suggested 
by WMAP measurements, early reionization should pro- 
duce a large-angle bump in the primordial GW-induced 
B-modes, allowing a possibly easier detection, without 
confusion with the CS-induced curl [19]. 

Actually, besides gravitational lensing, for cosmologi- 
cal models which constantly seed fluctuations in the ge- 
ometry, e.g. topological defects, vector metric perturba- 
tions can be huge and can prodiice non-negligible effects 
on the CMB photons as, in particular, B-mode polariza- 
tion, unlike to what happens in inflationary models [20]. 
On the other hand, no relevant contribution from these 
objects is indicated by the modern cosmological probes. 

In the present paper, we are interested in the cosmolog- 
ical stochastic GW background produced by Cold Dark 
Matter (CDM) halos via power transfer from scalar and 
possible vector perturbations to tensor metric modes, 
during the strongly non-linear stage of their evolution 
[21]. It differs from other cosmological backgrounds, as 
that produced during the mildly non-linear stage [22], 
since density and velocity fields can be, in this case, 
highly non-linear. 

Since the non-linear evolution of CDM halos occurs 
on a cosmological timescale, the produced gravitational 
radiation may be relevant at frequencies comparable 
to those of the primordial GW which affect the CMB 
photons and, therefore, can produce secondary CMB 
anisotropy and polarization, expecially B-modes, that 
could represent a foreground for the detection of the relic 
radiation. 

Moreover, as for the case of black holes and neutron 
stars, the analysis of the stochastic background produced 
by highly non-linear cosmic structures, could bring infor- 
mation on their distribution, evolution, shape and com- 
position, shedding light on many open issues. 

The plan of paper is as follows. In Sec. II we briefly out- 
line the mathematical formalism and show the analytical 
formulas we use to estimate the GW output from CDM 
halos. In Sec. Ill we introduce the homogeneous ellipsoid 
dynamics as an approximation to the halo virialization 
and adopt the halo mass function of Ref. [23] to describe 
their distribution. In Sec. IV we explain the technique 
used for the numerical evaluation of the stochastic GW 
background, while in Sec. V we show and discuss our 
results. Finally Sec. VI contains our concluding remarks. 



II. GRAVITATIONAL RADIATION: 
BASIC EQUATIONS 

The evolution of cosmological perturbations away from 
the linear regime is rich of several effects, such as the 



mode-mixing of different types of fluctuation, which not 
only implies that different Fourier modes influence each 
other, but also that density perturbations act as a source 
for curl vector modes and gravitational waves. 

Accordingly, cosmic structures can generate tensor 
metric modes during the non-linear stage of their evo- 
lution and, in particular, this mechanism applies to dark 
matter halos around galaxies and galaxy clusters in the 
highly non-linear regime. 

In the present paper, adopting the mathematical for- 
malism developed in Ref. [21], we estimate the output 
in gravitational waves from cosmic structures, following 
their evolution from the linear to the highly non-linear 
level. More specifically, the evaluation of this gravita- 
tional radiation is possible on scales much larger than 
the Schwarzscliild radius of collapsing bodies, by means 
of a "hybrid approximation" [21] of the Einstein field 
equations, which mixes post-Newtonian (PN) {e.g. [24- 
30]) and second-order perturbative techniques {e.g. [31- 
36]) to deal with the perturbations of matter and geom- 
etry. This approach gives a more accurate description of 
gravitational waves generated by non-linear CDM struc- 
tures than the standard second-order perturbation theory 
[22], which can only account for small deviations from 
the linear regime, or the Newtonian quadrupole radia- 
tion [37, 38]: indeed, it upgrades the weak- field limit of 
Einstein equations to account for PN scalar and vector 
metric, perturbations and for leading-order source terms 
of metric tensor modes. It provides, on small scales, a 
PN approximation to the source of gravitational radia- 
tion, and, on large scales, it converges to the first and 
second-order perturbative equations as obtained e.g. in 
Ref. [39], but still describing, on all the cosmologically 
relevant scales, the dynamics of the involved CDM struc- 
tures by means of the standard Newtonian Poisson, Euler 
and continuity equations {e.g. [40]) 

VV = 47rGa^^/9 , (1) 
p' + 3Wp + a,(pw")=0, (2) 

v'^ + HVa + V^d'^Va = -daf , (3) 

where is the gravitational potential associated with the 
density perturbation, p = p + Sp is the total matter den- 
sity composed by the background matter density, p, and 
the matter density perturbation, Sp, and, finally, v is 
the peculiar velocity field associated to the CDM halos. 
Greek indices denote spatial components; we adopt con- 
formal time r] and comoving coordinates x", in the Pois- 
son gauge, and assume that the Universe is spatially flat 
and filled with a cosmological constant A and a pressure- 
less fluid whose stress-energy tensor reads T^j = pu^Uj 
{u^Uj = — 1). Finally, H = a' /a, where primes indicate 
differentiation with respect to r] and a is the scale factor 
of the Universe which evolves according the Priedmann- 
Robertson- Walker background model. 
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Indeed, as the background cosmology, we have adopted 
a flat ACDM model with present baryon density given 
by fiob = 4.318 x 10-^ dark and CDM energy den- 
sity iloA = 0.7434, r^ocDM = 0.2134, Hubble constant 
Ho = 100ft km/sec/Mpc where h = 0.7199 and three 
massless neutrino species; the primordial perturbation 
spectrum is made by scalars only, normalized by cts = 
(j{R = 8 Mpc) = 0.9, with spectral index ris = 0.96 
[41, 42]. 

Accordingly to Ref. [21], in order to evaluate the 
stochastic background of gravitational radiation gener- 
ated by CDM halos, we will exploit the formula express- 
ing the solution of the inhomogeneous GW equation on 
scales well inside the Hubble horizon and in the so-called 
wave zone, which is 



(4) 



where r is the comoving distance between source and 
observer while the projection operator is given by V"' ^ = 
S'^fj — x°'xp/r'^. Eq. (4) expresses the GW output /?"^ 
in terms of integrals over the source "stress distribution" 
given by 



+ 4^ [d'^^dp^ - d,^ 6^^ . (5) 

The subscript "ref in Eq. (4) means that the quan- 
tity has to be evaluated at the retarded space-time point 
(// — r/c, x), i.e. at the source and at the emission time. 



III. THE ELLIPSOIDAL COLLAPSE MODEL 

Recently, N-body simulations in CDM models have 
shown departure of the halo density profile from the 
spherical symmetry {e.g. [43]) and suggest a triaxial 
shape which seems to be confirmed by optical. X-ray and 
lensing observations of galaxy clusters {e.g. [44, 45]). 

Consequently, according to the arguments in the pre- 
vious Section, CDM halos arc potential sources of gravi- 
tational radiation through power injection from the grav- 
itational potential and peculiar velocity, especially dur- 
ing the highly non-linear stage of their evolution, when 
density contrasts and velocity fields can be strongly non- 
linear. Since the aim of this paper is to evaluate the 
stochastic GW background generated by a distribution 
of cosmic structures, in this Section we will describe the 
model adopted to approximate their dynamics and viri- 
alization. 



A. The homogeneous elUpsoid dynamics 

We will use the gravitational collapse of homogeneous 
ellipsoids as described in Ref. [46], which developed a 



picture of cosmic structure formation that identifies viri- 
alized cosmological objects with peak patches in the ini- 
tial Lagrangian space. These peaks represent overden- 
sities in the initial Gaussian density field whose evolu- 
tion is approximated by a homogeneous ellipsoid dynam- 
ics. Each perturbation evolves under the infiuence of its 
own gravity and under the external tidal field (gener- 
ated by the surrounding matter) which, together with 
initial conditions, is chosen to reproduce the Zel'dovich 
approximation in the linear regime. Virialization is de- 
fined as the time when the third axis collapses and, fol- 
lowing Ref. [46], each axis is frozen once it has reached 
a freeze-out radius, chosen so that the density contrast 
at virialization, in the limit of spherical collapse, is the 
same as prescribed by the top-hat model. 

The peculiar velocity field is conveniently described in 
the system identified by the three principal axes, charac- 
terized by three different scale factors Ra {a = 1,2,3); 
thus, inside the homogeneous ellipsoid, peculiar velocities 
may be written as 



Ra . \ 
a— ajxa, 



(6) 



where we arc still adopting comoving coordinates but 
now time derivatives are with respect to the proper time 
dt = adr}. 

The internal peculiar gravitational potential, still with 
respect to the principal-axis system, is given by (see 
Ref. [46] for details) 



.a=l 



(7) 



where S = Sp/p is the matter density contrast while the 
factors ba are given by (see e.g. Refs. [47, 48]) 



R1R2R3 



{Rl + s),/{Ri + s) {Ri + s) {Ri + s) 

(8) 

Finally, the coefficients AwGa'^pX'^ are the eigenvalues 
of the traceless external tidal tensor (proportional to 
the traceless part of the peak strain) for which a lin- 
ear approximation is assumed [46], imposing that it 
evolves through the same equations satisfied by the linear 
growth-factor of density fluctuations in the considered 
cosmological background. 

After imposing the Zel'dovich approximation to fix the 
initial conditions on the proper ellipsoid axis lengths and 
their time derivatives, the evolution of an ellipsoidal per- 
turbation is specified through the equations [46] 



d^Rg _ hJ_ 
df^ ~ ~3" 



Ra - A-KGpRa 



+ 77 + ^0 + 



47rG _ Ac2 



(9) 



(10) 
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PR1R2R3 = const 



pa = const 



(11) 



(12) 



where, in Eq. (9), b'^ = - 2/3. 

Eqs. (9)-(12) have been numerically integrated using a 
fourth-order Runge-Kutta scheme and the integrals (8) 
have been evaluated by means of the so-called Carlson's 
elliptic function of the third kind. 

Fig. 1 shows the axis evolution versus time of a ho- 
mogeneous ellipsoid of mass A/ = 5 x lO^^AfQ and initial 
overdensity 6{zi = 40) = 6.4 x 10~^ at comoving distance 
D = 100 Mpc from the observer. The shape of the el- 
lipsoid is the most probable in terms of the distribution 
of ellipticity and prolateness, to be defined in the next 
sub-Section. The evolution follows Eqs. (9-12) and the 
axis freezing out method suggested in Ref. [46]. Let us 
stress that, contrary to other ellipsoidal collapse schemes 
{e.g. [49]), this model implies that virialization is reached 
when the third and not the first axis collapses, while the 
freezing out method avoids d — > 00. 

In order to estimate the GW output by CDM cosmic 
structures, we insert Eqs. (6)-(7) in Eqs. (4)-(5) during 
the collapse of each homogeneous ellipsoid which repre- 
sents, in our simulation, a CDM halo evolving towards 
virialization. 

In Fig. 2 we show two of the three non- vanishing trace- 
less source components generated by the halo collapse 
of Fig. 1. These components are evaluated with re- 
spect to the eigenframe of the ellipsoid principal axes 
at rest with respect to the expanding cosmological back- 
ground; by performing a transverse projection, the grav- 
itational waves in the observer frame are obtained. Ac- 
tually, Fig. 2 represents these two components divided 
by 1 -I- 2; = l/a(i), in order to separate the effects of the 
background expansion, included in Eq. (4) , from the halo 
evolution itself. 



B. The most probable ellipsoid and the halo mass 
function 



Once the cosmological background model is fixed, the 
evolution of an ellipsoidal perturbation is determined by 
three parameters given by the three initial eigenvalues 
of what, in the Zel'dovich approximation, is called the 
deformation tensor, da/s — (l/a^)V^V^(/3; the latter are 
related to the initial ellipticity e, prolateness p and linear 
density contrast S of the perturbation; those relations 
read [23] 



P = 



Ai — A3 
26 ' 

Ai + A3 — 2A2 



Ai -f A2 4- A3 



(13) 



(14) 



(15) 
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FIG. 1: Evolution of the principal axis scale factors for the 
most probable ellipsoid of mass M = 5 X 10^^ M© and ini- 
tial overdensity 5{zi = 40) — 6.4 x 10~^, embedded in a flat 
ACDM universe, at distance D = 100 Mpc from the observer. 



where the Aq are the eigenvalues of dap with Ai > A2 > 
A3, which, if (5 > 0, implies e > and — e < p < e. 

For a Gaussian random density field, smoothed in real 
space with a top- hat filter of size V = 47ri?'^/3 and mass 
M — AnpR^ /?>, on average and for a given 5 the prolate- 
ness is p = 0; consequently, the most probable ellipticity 
is e„jp = {a/S)/^/5. Here a ~ cr{R) represents the linear 
rms value of the S distribution [23]. 

From these considerations and from the homogeneous 
ellipsoid collapse model as described in Ref. [46] , the au- 
thors of Ref. [23] have determined the shape of the mov- 
ing barrier, i.e. the critical overdensity required for CDM 
structure virialization at redshift z: that is 



(16) 



where v = [5sc{z)/a{M)] , Ssdz) is the critical overden- 
sity required for spherical collapse at z extrapolated using 
linear theory to the present time, and cr is the linear rms 
value of the initial density fluctuation field also extrapo- 
lated to the present time. The parameters (3 « 0.485 and 
a w 0.615 come from ellipsoidal dynamics and the value 
q ~ 0.75 comes from normalizing the model to simula- 
tions [50]. 

Using Eq. (16) in the excursion set approach in order to 
obtain the distribution of the first crossings of the barrier 
by independent random walks, the authors of Refs. [23, 
51] have derived the average comoving number density 
of halos of mass M, i.e. the so-called unconditional halo 
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where n 
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FIG. 2: Two of the three non- vanishing traceless source com- 
ponents generated by the halo collapse of Fig. 1 



mass function 
n(M, z)dM = 



■K M2 cr(M) 
diner 

exp 



2p 



dlnM 



2a{Mf 



dM , (17) 



where po is the mean comoving cosmological mass den- 
sity, while p = 0.3 and A = 0.32218. The Press-Schechter 
mass function is recovered for q = 1, p — and A — 0.5 
[52]. 

In what follows, a{M) and Ssc{z) are computed ac- 
cording to the formulas [53-55] 

0- cx (1 + 2.208m'^ - 0.7668m^'^ + 0.794W'^)-^/(^'^) , 

(18) 

where d = 0.0873, m = M{rh)^ and 

r = r!o,„/iexp[~riob(i + \/2^/rJom)] . (i9) 

The quantities related to the density contrast are 

S, D+iz = 0) 



D+{z) 



, 3(127r)2/3^ 

"~ 20 (1 + 0-0123 logio an) 



(20) 



(21) 



The linear growth factor of density fluctuations, normal- 
ized to unity at present, may be approximated as [56] 



2(1 + z) 



On 



a„r-l7A+ 1 + ^ 1 + i^ 



70 



1 + zf/E'^{z), f^A = ^gk/E'^{z), and 
E{z)^H{z)/Ho = [nora{l + zf 



OA 



,1/2 



(23) 



IV. THE STOCHASTIC GW BACKGROUND 

In order to evaluate the GW output generated by a 
spatial distribution of CDM halos we will exploit Eq. (17) 
which provides a good fit to N-body simulations of struc- 
ture clustering in a variety of cosmological models, at 
least over the redshift range z = — 4 [23, 50, 51, 57, 58]. 

For theoretical consistency, we have chosen to follow 
the same strategy adopted by Ref. [23] as described in 
the previous section. Therefore, in our numerical com- 
putation, we consider CDM structures over a mass range 
Af = 5 X IO^Mq - 5 X IO^^Mq, which viriahze at redshifts 
from z = to 4. Each of these structures is approximated 
by a homogeneous ellipsoidal perturbation with mass A/, 
linear mass variance (T^(Af) and critical linear density 
contrast 5{M,z) = B{a^,z); in other words, every per- 
turbation represents the most probable ellipsoid (p = 
and e — Cmp) of mass M which collapses at redshift z. 

Given the density contrast, the ellipticity and the pro- 
lateness, we then calculate the eigenvalues of the ex- 
ternal tidal tensor, using Eqs.(13)-(15) and the relation 
= Aq — (5/3. Next, we linearly rescale all quantities 
to the initial redshift Zi = 40, at which the ellipsoidal 
evolution of the density perturbation starts, following 
Eqs. (9)-(12). In fact, while the mass function provides 
the number of halos virializing at a given redshift (in our 
case z = — 4) , the evolution of matter density perturba- 
tions, giving rise to these virialized objects, begins much 
before, i.e. at very high redshifts (in our case Zi = 40). 
The initial conditions on the scale factor are given by the 
relation a{zi) = 1/(1 -I- Zi) and by the well-known Fried- 
mann equations, while, as we have already anticipated, 
the initial conditions on the axis lengths and their time 
derivatives are specified by the Zel'dovich approximation 
setup 

R^{zi) ^ a{z,){l ~ X^) (24) 



and 



i?„(z,) - Hiz,) [i?„(z,) - a(z,)/(z,)A„] , (25) 



(22) 



where f{z) « fl^f + (1/70) [1 - l/2f7m (1 + f^m)] is the 
growth rate of density fiuctuations {e.g. Ref. [59]). 

For each M and z, using Eq. (4) and switching to the 
proper time t, we evaluate the two independent compo- 
nents of the gravitational radiation produced by a CDM 
halo, assuming that it is casually oriented and placed at a 
comoving distance r(z) from the observer, where z is the 
collapse redshift. In this way, we observe today the ra- 
diation emitted at the virialization time when, according 
to our ellipsoidal model, the GW output has the maxi- 
mum value. Actually, adopting this strategy, we slightly 
underestimate the total GW background, since we do not 
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take into account those CDM halos which are still away 
from virialization. Moreover it is noteworthy here that, 
in our approach, we have extrapolated Eq. (4) outside its 
range of validity. In fact, this formula holds on scales well 
inside the Hubble horizon and in the wave zone {i.e. at 
distances larger than both the characteristic wavelengths 
and the characteristic size of the source), while, as we pre- 
viously noticed, CDM halos generate gravitational waves 
whose frequency is comparable with the inverse of the 
Hubble time. Nonetheless, as shown in the next Section, 
our results agree with several analytic approximations 
and previous works. 

To account for all the directions of observation, we con- 
vert the two independent states /i"^ (t) of tensor polariza- 
tion from the frame associated with the ellipsoid principal 
axes to the observer frame, assuming that CDM struc- 
tures emit in all directions and are uniformly distributed 
all around the observer. For this purpose, we use the rela- 
tion h°'i}{t,n) = i?'^«^(f])/i^^(i)i?'^^(n) where i?"/3(f^) 
is the general form of the rotation matrix with = [60] 

cos sin ^ 

R°'/s{0,)= I — cos0sin0 cos6'cos(/) sin 6' 
sin 6 sin (f> — sin 9 cos (f> cos 9 

and the solid angle f2 = {9,(j)) is defined following the 
conventions of Ref. [61]. 

Since our aim is to estimate the energy density 



^GW(J^) 



3Hi 



(26) 



associated with the stochastic GW background at the ob- 
server {e.g. Ref. [1]), we need to know the power spectral 
density PSD(:/), which one can obtain from the Parseval's 
theorem as 



di^FSB{iy) . 



(27) 



That depends on the redshifted proper frequency u = 
Vf,/{l+z), where is the proper frequency at the emis- 
sion time. In Eq. (27) angle brackets denote time aver- 
aging at a given spatial point. 

Thus, we first numerically evaluate the 
Y'ST>{v,z,M,Q) of each individual component of 
h°'f3{t,fl) at each fixed value of M, z and fl, then 
we average the calculated PSDs over all directions by 
integrating over the solid angle and dividing by 47r and, 
finally, we sum over the components in order to get a 
mean power spectral density PSD{i^, z, M) for every z 
and M. 

Since in our model each CDM halo is approximated by 
a most probable ellipsoid of mass M which collapses at 
redshift z, we multiply each PSD(i^, z, M) by the number 
dN{z, M) = n{z, M)dMdV of halos in the comoving vol- 
ume dV{z) = 4:C'Kr^{z)dz/ {aoH{z)) where c is the speed 
of light and ao = 1 is the present value of the scale factor. 
Finally, we insert the resulting quantity in the definition 



of the GW energy density Eq. (26) and integrate over all 
redshifts and masses to obtain the total f2Gw(i^)- 
All the results are presented in the next section. 



V. RESULTS 

Our result concerning the GW output of each CDM 
halo, an example of which is given in Fig. 2, is consistent 
with previous works in this field [62, 63]. Moreover it is 
comparable to analytic approximations {e.g. [37, 64]) as 



3 X 10-" GM (10^5) 
D/lOOMpc ^ 



2/7 



M 



1O15M0 ' 



(28) 



where h represents the amplitude of a GW signal coming 
from a non-spherically symmetric collapsing object with 
characteristic size L at distance D from the observer. 

It is worth noting that the produced gravitational radi- 
ation has a very long characteristic period, approximately 
given by the inverse of the halo evolution time, which, ac- 
cording to the ellipsoidal model, corresponds to frequen- 
cies of the order of « 10-^^Hz. This excludes, there- 
fore, any direct detection of a complete pulse, but still 
allows for the possibility of GW detection via secondary 
CMB anisotropy and polarization and via the "secular 
efii'ect" discussed in Refs. [62, 63]. The latter takes place 
when a gravitational-wave crosses two testing particles; 
this induces a variation in their relative distance which 
increases in time, since this efi'ect lasts for many years. 

Actually, besides what stressed in the previous Sec- 
tion, there are other reasons for which the ellipsoidal 
collapse approximation to CDM halo virialization un- 
derestimates amplitude and frequency characterizing the 
GW background. In fact, using this approach, the evo- 
lution of cosmic structures is regarded as a continuous 
phenomenon which neglects merging effects and any pos- 
sible features of variability that, according to Ref. [62], 
should be characterized by a dynamical frequency of the 
order of w IQ-'^^Hz. 

In Fig. 3 the main result of our paper is shown, i.e. the 
total energy density r2Gw('^) ~ 10"^", associated with 
the stochastic halo-induced GW background, as a func- 
tion of the proper frequency u at the observation. The 
total spectrum of the signal is composed by many single 
peaks which represent the contribution to the total back- 
ground from each most probable halo weighted via the 
mass function at different redshifts. On the other hand, 
as the following discussion shows, these peaks are caused 
by the subset of structures leading to a non-negligible 
GW signal. In fact Eq. (4) shows that the GW amplitude 
is proportional to the inverse of the comoving distance, 
while, from the expression of the eSiciency e = GM/{c^L) 
and the total radiated energy Eqw = eMc^, where L rep- 
resents the characteristic halo size at virialization, it fol- 
lows that more massive objects give rise to higher values 
of the GW strain. This effect is also confirmed by numer- 
ical estimates of the power spectral density for different 
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FIG. 3: The total energy density Qowii') associated with 
the stochastic GW background induced by CDM halos as a 
function of the proper frequency v at observation. 




FIG. 4: Contribution to the total energy density fiGwC^^) by 
two most probable CDM halos of mass 5 x 10^^ A/© placed at 
redshifts z = 0.025 and 0.075. 



objects in the redshift range < z < 4. In fact, masses of 
the order of 10® — lO^M©, although weighted via the mass 
function in Eq. (17), contribute to Qqw{^) by only a fac- 
tor of orders of 10"'^'^ — 10"^*; since the amplitude of the 
gravitational waves decreases with distance, the greater is 
the redshift z, the lower is their contribution. Thus only 
a few peaks are visible in Fig. 3 since the energy density 
produced by less massive structures is completely negligi- 
ble with respect to the effect (of orders of 10~^^ — 10^^") 
of far more massive objects (10" - lO^^M©) at low red- 
shifts, z < 1. Consequently, the dominant contribution 
to the stochastic GW background is likely to be pro- 
duced by CDM halos corresponding to nearby galaxies 
and galaxy clusters which contribute by several orders 
of magnitudes more than their substructures, although 
the latter are far more numerous. Indeed, in Fig. 4 we 
may look at the contribution to the total Ogvv('^) (see 
Fig. 3) by two most probable halos of mass 5 x 1O^^M0, 
placed at redshifts z — 0.025 and 0.075. In its maximum 
height, the signal reaches about half of the correspond- 
ing value in Fig. 3. The remaining part of the signal is 
caused by many halos of comparable mass, as well as by 
those about one order of magnitude lighter, for which 
the mass decrease is compensated by the increase in the 
number. It is worth noting that the ellipsoidal collapse 
model introduces a three-peak pattern due to the freez- 
ing out method used to stop the axis collapse, which as 
a zero-th order approximation imposes stability at viri- 
alization, ignoring any residual dynamics. In the case of 
the specific geometrical configuration of the most prob- 



able ellipsoid considered in Fig. 4, this translates in two 
prominent peaks and a third negligible one. Actually, 
the residual dynamics at virialization would most prob- 
ably imply a broadening of the spikes, decreasing the 
frequency splitting, possibly converging to a single peak 
for configurations close to sphericity. 

Finally and most importantly, the quantity 
h^nQ\fli{i/) w 10"^" is comparable to the energy 
density associated with the stochastic background 
induced by primordial GWs. In fact, if the energy scale 
of inflation is V^^"^ « 1 - 2 x lO^^CeV, the energy 
density associated with the primordial stochastic GW 
background, with a tensor spectral index nx ~ 0, 
is « 10^^^ — lO"""^^ for frequencies of the order of 
10-1* - 10"" Hz (e.g Ref. [1]). 



VI. CONCLUDING REMARKS 

In this work, we have estimated the GW background 
from cosmological tensor modes produced by the highly 
non-linear collapse of CDM density perturbations, i.e. 
generated during the strongly non-linear stage of CDM 
halo evolution. 

We found that the signal is significant at very low fre- 
quencies, « 10^1* Hz, as a consequence of the cos- 
mological time scales involved in the collapse of CDM 
halos. This signal appears as a broad peak made by the 
superposition of many impulses, all centered around fre- 
quencies of the order of 10-^* Hz. Most importantly, our 
results suggest that the signal is likely to be comparable 
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to the primordial tensor power if inflation occurred at the 
GUT scale. 

We want to stress that the homogeneous ellipsoidal 
collapse model, adopted to simulate CDM halo evolution 
and virialization, underestimates the frequency and am- 
plitude of the emitted gravitational waves, since, at each 
redshift z, it does not take into account non-virialized 
objects and neglects variability features and merging ef- 
fects that could enhance the anisotropic stress sourcing 
tensor modes, which are more sensible to the velocity field 
rather than to the peculiar gravitational potential. Con- 
sequently, the total energy density riQ^iy) generated by 
cosmic structures could even be of one or two order of 
magnitudes greater and overcome the stochastic back- 
ground associated with primordial gravitational waves at 
the same frequencies (see also results in Ref. [63]). 

The CDM halo GW background could also produce a 
non-negligible contribution when considering the cosmo- 
logical tensor-to-scalar ratio. 

Due to the cosmological scales involved, and to the 
amplitude of the signal, it is reasonable to expect that 
these gravitational waves could affect the primary CMB 
anisotropics, contributing to the Integrated Sachs Wolfe 
(ISW) effect caused by the time evolution of cosmologi- 
cal perturbations between us and the last scattering sur- 
face. The stochastic GW background from CDM halos 
might boost the temperature anisotropics on large angu- 
lar scales, where however the contribution from density 
fluctuations dominates. On the other hand, the produced 



temperature quadrupole can be scattered off by the free 
electrons of the intra-cluster and intra-galactic media, 
giving rise to secondary E and B polarization modes sim- 
ilarly to what happens for the primordial temperature 
quadrupole as described in Ref. [65] . These contributions 
have to be taken into account when performing a precise 
evaluation of the level of CMB polarization anisotropy 
expected for the forthcoming polarization oriented CMB 
probes, in particular for what concerns the B modes. 

Most of these issues deserve a careful investigation in 
future works. Here we conclude stressing again our main 
results, suggesting that the amplitude of the stochastic 
GW background generated by CDM halos in their non- 
linear evolutionary phase is comparable or larger than the 
signal expected from the early universe in the inflationary 
scenario. We also remark that our findings are consistent 
with existing analytical approximations. The forthcom- 
ing steps are the improvement of the calculation of the 
source of the signal, making use of cosmological N-body 
simulations, as well as the computation of the induced 
CMB anisotropy in total intensity and polarization. 



Acknowledgments 

We warmly thank Bepi Tormen for helpful discussions 
and suggestions. 



[1] M. Maggiore, Phys. Rep. 331, 283 (2000). 

[2] V. Ferrari, S. Matarrese, R. Schneider, Mon. Not. Roy. 

Astron. See. 303, 258 (1999). 
[3] K. A. Postriov, in Free. "Particles and Cosmology", 9th 

International Baksan School, 1997, [astro-ph/9706053]. 
[4] D. I. Kosenko and K. A. Postnov, Astron. Astrophys. 

336, 786 (1998). 
[5] R. Schneider, A. Ferrara, V. Ferrari, S. Matarrese, Mon. 

Not. Roy. Astron. Soc. 317, 385 (2000). 
[6] R. Schneider, V. Ferrari, S. Matarrese, S. F. Porte- 

gies Zwart, Mon. Not. Roy. Astron. Soc. 324, 797 (2001). 
[7] A. J. Farmer and E. S. Phinney, Mon. Not. Roy. Astron. 

Soc. 346, 1197 (2003). 
[8] L. P. Grishchuk, [arXiv:gr-qc/0305051]. 
[9] W. HU, U. Seljak, M. White, M. Zaldarriaga, Phys. Rev. 

D 57, 3290 (1998). 
[10] M. Kaniionkowski and A. Kosowsky, Ann. Rev. Nucl. 

Part. Sci. 49. 77 (1999); W. Hu and S. Dodelson, Ann. 

Rev. Astron. Astrophys. 40, 171 (2002). 
[11] S. Dodelson, Modem Cosmology (Academic Press, 2003); 

A. Liddle and D. H. Lyth, Cosmological Inflation 

and Large-Scale Structure (Cambridge University Press, 

2000). 

[12] J. Pritchard and M. Kamionkowski, Annals Phys. 318, 
2 (2005). 

[13] M. Kamionkowski, A. Kosowsky, A. Stebbins, Phys. Rev. 

D 55, 7368 (1997). 
[14] M. Zaldarriaga and U. Seljak, Phys. Rev. D 55, 1830 



(1997). 

[15] U. Seljak and M. Zaldarriaga, Phys. Rev. Lett. 78, 2054 

(1997) . 

[16] M. Bartelmann and P. Schneider, Phys. Rept. 340, 291 
(2001). 

[17] M. Zaldarriaga and U. Seljak, Phys. Rev. D 58, 023003 

(1998) . 

[18] U. Seljak and C. M. Hirata, Phys. Rev. D 69, 043005 

(2004) . 

[19] P. Cabella and M. Kaniionkowski, Lectures given at 
the 2003 Villa Mondragone School of Gravitation and 
Cosmology: "The Polarization of the Cosmic Mi- 
crowave Background," Rome, Italy, September 6-11, 
2003, [arXiv:astro-ph/0403392]. 

[20] N. Turok, U. L. Pen, U. Seljak, Phys. Rev. D 58, 023506 
(1998). 

[21] C. Carbone and S. Matarrese, Phys. Rev. D 71, 043508-1 

(2005) . 

[22] S. Matarrese and S. Mollerach, "The stochastic 
gravitational-wave background produced by non-linear 
cosmological perturbations", in Proc. "Some topics on 
General Relativity and Gravitational Radiation", p. 87, 
edited by J. A. Miralles, J. A. Morales, D. Saez, Editions 
Frontires, Paris, 1997, [arXiv:astro-ph/9705168]. 

[23] R. K. Sheth, H. J. Mo, G. Tormen, Mon. Not. Roy. As- 
tron. Soc. 323, 1 (2001). 

[24] S. Chandrasekhar, Astrophys. J. 142, 1488 (1965). 

[25] S. Chandrasekhar, Astrophys. J. 158, 45 (1969). 



9 



S. Chandrasokhar, Astrophys. J. 160, 153 (1970). [49 
K. Tomita, Prog. Thcor. Phys. 79, 258 (1988). [50^ 
K. Tomita, Prog. Thcor. Phys. 85, 1041 (1991). 
M. Shibata, H. Asada, Prog. Theor. Phys. 94, 11 (1995). [51 
S. Mataxrese and D. Terranova, Men. Not. Roy. Astron. 
Soc. 283, 400 (1996). [52 
T. Pyne and S. Carroll, Phys. Rev. D 53, 2920 (1996). 
S. MoUerach and S. Matarrese, Phys. Rev. D 56, 4494 [53 
(1997). [54 
S. Matarrese, S. MoUerach and M. Bruni, Phys. Rev. D 
58, 043504 (1998). [55 
K. Tomita, Prog. Theor. Phys. 37, 831 (1967). 
S. Matarrese, O. Pantano and D. Saez, Phys. Rev. Lett. [56 
72, 320 (1994). 

K. Tomita, Prog. Theor. Phys. 45, 1747 (1971); 47, 416 [57 
(1972). 

C. W. Misner, K. S. Thorne and j. A. Wheeler, Gravita- 
tion, W. H. Freeman and Company, San Francisco (1973). [58 
L. D. Landau and E. M. Lifshitz, The classical theory of 
fields, Pergamon Press, New York, (1975). 

S. MoUerach, D. Harari and S. Matarrese, Phys. Rev. D 
69, 063002 (2004). [59 
P. J. E. Peebles, The Large-Scale Structure of the Uni- 
verse, Princeton University Press, Princeton (1980). [60 
M. Tegmark et al, Phys. Rev. D 69, 103501 (2004). 

D. N. Spergel et al., Astrophys. J. Suppl. 148, 175 (2003). [61 
Y. P. Jing and Y. Suto, Astrophys. J. 529, L69 (2000). [62 
M. J. West, Astrophys. J. 347, 610 (1989). 

M. Plionis, J. D. Barrow, C. S. Frenk, Mon. Not. Roy. [63 
Astron. Soc. 249, 662 (1991). 

J. R. Bond and S. T. Myers, Astrophys. J. Suppl. Series [64 
103, 1 (1996). 

S. Chandrasekhar, Ellipsoidal Figures Equilibrium, Yale 
Univ. Press, New Haven, (1969). [65 
J. Binney and S. Tremaine, Galactic Dynamics, Prince- 
ton Univ. Press, Princeton, (1987). 



S. D. White and J. Silk, Astrophys. ,L 231, 1 (1979). 
R. K. Sheth and G. Tormen, Mon. Not. Roy. Astron. Soc. 
329, 61 (2002). 

R. K. Sheth and G. Tormen, Mon. Not. Roy. Astron. Soc. 
308, 119 (1999). 

J. R. Bond, S. Cole, G. Efstathiou, N. Kaiser, Astrophys. 
J. 379, 440 (1991). 

T. Kitayama and Y. Suto, Astrophys. J. 469, 480 (1996). 
J. A. Peacock and S. J. Dodds, Mon. Not. Roy. Astron. 
Soc. 267, 1020 (1994). 

N. Sugiyiama, Astrophys. J. Suppl. Series 100, 281 
(1995). 

S. M. Carroll, W. H. Press, E. L. Turner, Annu. Rev. 
Astron. Astrophys. 30, 499 (1992). 

G. Kauffmann, M. J. Colberg, A. Diaferio, 
S. D. M. White, Mon. Not. Roy. Astron. Soc. 303, 
188 (1999). 

A. Jenkins, C. S. Frenk, F. R. Pearce, P. A. Thomas, 
M. J. Colberg, S. D. M. White, H. M. P. Couchman, 
J. A. Peacock, G. Efstathiou, A. H. Nelson, Astrophys. 
J. 499, 20 (1998). 

O. Lahav, P. B. Lilje, J. R. Primack, M. Rees, Mon. Not. 
Roy. Astron. Soc. 251, 128 (1991). 

H. Goldstein, Classical Mechanics, Addison- Wesley, 
Cambridge, Mass., (1956). 

R. L. Forward, Phys. Rev. D 17, 379 (1978). 

V. Quilis, J. M. Ibanez and D. Saez, Astrophys. J. 501, 

L21 (1998). 

V. Quilis, J. M. Ibanez and D. Saez, Astron. Astrophys. 
353, 435 (2000). 

S. L. Shapiro and S. A. Teukolsky, Black Holes, White 
Dwarfs ans Neutron Stars, John Wiley, New York, 
(1983). 

G. C. Liu, A. da Silva, N. Aghanim, Astrophys. J. 621, 
15 (2005). 



